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Abstract: We derive the Leading-Order master equation to extract the polarized gluon 

^ I distribution G{x,Q'^) = x5g{x,Q'^) from polarized proton structure function, g\{x.,Q'^). 

By using a Laplace-transform technique, we solve the master equation and derive the 



f"^ ■ polarized gluon distribution inside the proton. The test of accuracy which are based on 



our calculations with two different method confirms that we achieve to the correct solution 
for the polarized gluon distribution. We show that accurate experimental knowledge of 

^^ \ g^{x , Q"^) in a. region of Bjorken x and Q^, is all that is needed to determine the polarized 

gluon distribution in that region. Therefore, to determine the gluon polarization -^ ,we 
only need to have accurate experimental data on un-polarized and polarized structure 

X ■ functions {F^{x, Q^) and c/f (x, Q^)). 
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1. Introduction 

The spin structure of the proton is one of the most challenging open puzzles in Quantum 
Chromodynamics. One of the main questions in particle and nuclear physics is: How is the 
proton spin built up from its quark and gluon constituents? Now, we know that the quark 
contribution to the spin of the proton is about 0.30. It is determined precisely in a QCD fit 
to the polarized proton structure function data, g\{x,Q'^). Therefore at the present time, 
the role played by the gluons in the nucleon spin is the most challenging task. 
The main goal of the experiments at HERMES, SMC, COMPASS and RHIC spin physics 
programs is the measurement of the helicity contribution of the gluons to the nucleon spin, 
Ag. Experimentally this value is mainly accessible via two processes in polarized DIS ex- 
periments. The first one is the production of high pt hadron pairs with large transverse 
momentum, but these processes have large background contributions from QCD Compton 
processes and fragmentation that one has to control. The second one is open charmed 
meson production in photon-gluon fusion process. The cross section of these two processes 
is directly related to the ratio of the polarized gluon density to the un-polarized gluon 
density, (q'1\ [Hi 0]- Then, to measure the gluon contribution to the nucleon spin with 
this approach, it is necessary to know the un-polarized gluon distribution well. Polarized 
gluon measurements from deep inelastic experiments are summarized in Table 1. 
In line with experiments, many NLO QCD global fits to the inclusive ^^(x, Q^) data were 
used to extract the magnitude of /S.g and the shape of 6g{x,Q'^) |^, ^, |lO|. Unfortunately, 



Experiment 


process 


<Xg> 


< fi"^ > 




Reference 


HERMES 


hadron pairs 


0.17 


~2 


0.41 ±0.18 ±0.03 


i 


HERMES 


inclusive hadrons 


0.22 


1.35 


0.071 ±0.034 ±0.105 


i 


SMC 


hadron pairs 


0.07 


— 


-0.20 ±0.28 ±0.10 


i 


COMPASS 


hadron pairs, Q^ < 1 


0.085 


~3 


0.016 ±0.058 ±0.055 


ii 


COMPASS 


hadron pairs, Q^ > 1 


0.13 


— 


0.06 ±0.31 ±0.06 


ii 


COMPASS 


open charm 


0.15 


13 


-0.57 ±0.41 ±0.17 


i 



Table 1: Polarized gluon measurements from deep inelastic experiments. 



these global fits are sensitive to the initial assumption for the polarized gluon distribution 
and also to the initial scale of Qq. The shape of the polarized gluon distribution extracted 
from different groups are not identical because of their different assumptions for the initial 
polarized gluon distribution, initial Qq and also different approaches for global fits. We 
would like to know how big the gluon spin contribution is to the total spin of the proton 
and look for a simple way to measure this value. This work can be considered as a proposal 
for direct measurement of the polarized gluon distribution inside the proton with more ac- 
curacy. 

Recently an explicit expression for the un-polarized gluon distribution function G{x, Q^) = 
xg{x,Q'^) in the proton in terms of the proton structure function ^^'(a;,^^) was derived 

Here, the same procedure is used to 
xSg{x,Q^), inside the pro- 



by using a Laplace-transform technique |11, 12, |i| 
derive the polarized gluon distribution function, G{x,Q'^) 
ton. We obtain an analytic solution for the polarized gluon distribution in terms of the 
polarized proton structure function, g^{x,Q'^). Thus we can calculate the polarized gluon 
distribution inside the proton directly only by finding g^^x, Q'^) over large kinematic range 
of X and Q^. This is the principal theoretical result in this paper. 

This paper is organized as follows. In section 2, we derive LO master equation for extract- 
ing the polarized gluon distribution G(x, Q"^) = x5g{x, Q^) from polarized proton structure 
function, g\{x,Q'^). Section 3 is devoted to the solution of the master equation. In the 
last section, we first perform a global parametrization of gi{x,Q'^) using all available ex- 
perimental data and then we calculate numerically the polarized gluon distribution. Our 
conclusion is given in section 5. 



2. LO master equation to extract the polarized gluon distribution G{x, Q^) 

x5g{x^Q'^), using the polarized proton structure function, g\{x,Q'^) 



The LO DGLAP equation, [14, 15, |T^ for the evolution of the polarized proton structure 
function, g\{x,Q'^) can be written as (See Appendix A) 



dglix^Q"^) as [^ dz p ^2-,,r^ (^ 



as_ 

'2^ 



E' 



''^'g{z,Q^)6K, '"^ 



^2^,.,^ ,„..,g^^ 



(2.1) 



where 5Kqq{x) and 5Kqg{x) are the LO polarized sphtting functions and Us is the renor- 
mahzed running coupling constant. We introduce Qi{x,Q^) as 

e.(,.Q=, = ,Mi£|!) _ ^,£ ^,r,,g=,«- , (£) , ,2.2, 

Now we can write the master equation, the DGLAP equation for the g^{x,Q'^), as follows 

gg{x, Q2) = X r ^Giz, Q^)5Kqq{-), (2.3) 

where gg{x,Q'^) = {^T,e^i)-^g^{x,Q^). In the Eq. (2.2) and Eq. (2.3), the LO q -> q 
and g ^ q polarized splitting function are given by [|ll 



6Kqqix) = ^((\^)+ + l^ix - 1)), (2.4) 



1. ^ 1 

2(2x-l) = x--. 



5/s:gg(x) = -(2x-l) = x--. (2.5) 



To extract the polarized gluon distribution, G{x, Q'^) = x6g{x, Q'^) from the DGLAP equa- 
tion, we should solve the master equation (2.3) and find the polarized gluon distribution 
inside the proton. This issue is the subject of the next section. 

3. Solution of the master equation 

To solve the master equation (2.3), we follow the procedure that was used by BDM |11| 
and used the Laplace transformation to solve this equation. Now we use the coordinate 
transformation as 

V = ln(l/x), (3.1) 

Then the functions G, Kqg, and Q in u-space are given by 

G{v,Q')^Gie'^,Q'), (3.2) 

dKqq{v,Q^) = 6Kqq{e-\Q^), (3.3) 

Giv,Q^)^gg{e-'',Q^). (3.4) 

Explicitly from Eq. (2.5), we have 

5Kqg{v,Q^) = e-''-^. (3.5) 

Therefore 



g{v,Q^)= G{w,Q^)e-'-''-'"^6kqg{v-w)dw= G{w,Q^)H{v - w)dw, (3.6) 
Jq Jo 



where w = ln{l/z) and H{v) is defined as 



H{v) ^ e-'"6K,g{v) = e-2- - ^e'^ (3.7) 

If we take the laplace transform of Eq. (3.6), then we have 

C[g{v,Q^);s]=C[[ G{w,Q^)H{v-w)dw;s] (3.8) 

Jo 

=^ Qis, Q2) = Gis, Q2) X his). (3.9) 

In the above equation we used the following property for Laplace transformation 

C~\F{s)G{s)] = f fit- T)giT)dT = [ fiT)git - T)dT. (3.10) 

Jo Jo 

So we have the polarized gluon distribution in s-space as 

Gis,Q') = ^^^. (3.11) 

Now the polarized gluon distribution in v-space is given by 

Giv,Q^) = C~^[gis,Q^)his)~^]= I giw,Q^)Jiv-w)dw, (3.12) 

Jo 

Where Jiv) = C^^[his)^^,v]. The calculation of J(f) by using the Eq. (3.7) and 
inverse Laplace transform of /i(s)~^ for LO is straightforward and given in term of the 
Dirac delta function 

Jiv) = 4 + 66iv)+26'iv). (3.13) 

Using the Eq. (3.12) and Eq. (3.13) the polarized gluon distribution in v-space is 
given by 

pv rv pv 

Giv,Q^) = i giw,Q^)dw + 6 5iv-w)giw,Q^)dw + 2 6' iv - w)giw,Q^)dw. 
Jo Jo Jo 

(3.14) 

Using the following relation for Dirac delta function 

rv 

/ 5iv-w)fiw)dw = fiv), (3.15) 

Jo 

pv 

/ 5'iv-w)fiw)dw = f'iv). (3.16) 

Jo 

we will have 



Giv, q2) = 4 r giw, Q^)dw + egiv, q^) + 2^^!^. (3.17) 

Jo dv 



Finally we have the polarized gluon distribution as 



Gix, Q2) = 4 ^' gg{z, Q^)^ + 6gg{x, Q^) - 2(x^^^^^), (3.18) 

where gg{x,Q'^) = {^T.4)~^Qi{x^Q'^) and Ql{x,Q'^) is given by Eq. (2.2). We use 4 
massless quarks (u,d,s,c) in our calculation and then '^e^ = ^- 

To calculate the right hand of the Eq. (3.18), we have to have an analytic function for 
g\{x^ Q^). So in the next section we try to do the global parametrization of g\{x, Q^) using 
all available experimental data and find this analytic function and then we calculate the 
polarized gluon distribution function numerically. 

4. Numerical results 

In this section we intend to use Eq. (3.18) to extract the polarized gluon distribution inside 
the proton. 

4.1 Global parametrization oi g\{x,Q'^) using all available experimental data 

We have parameterized the polarized proton structure function, g^{x,Q^) in < x < 1 as 

xglix, Q^) = x^{l - x) \ y\ J l + A{Q')Ln[^^j^^]). (4.1) 

\^ Xp) Xp X\i. Xp) 

Here xp = 0.234 is an approximate fixed point observed in the data where the curves for 
different Q^ cross. At that point, dg^{xp,Q'^)/dlnQ'^ ^ for all Q^; gp = g{{xp,Q'^) = 
0.241 is the common value of g^. We used all available experimental data for polarized 



proton structure function from E143, SMC, HERMES 2006 and COMPASS 2009 |T|, |19 



|20| , 21, |2^. More accurate global fit with more data, over large kinematic ranges of x and 
Q^ will lead to the precisely determination of the polarized gluon distribution function. 
The Q^ dependence of gi{x,Q'^) in our global fit is given by 

A{Q'^) = ao + ai\nQ'^ + a2\n^Q'^. (4.2) 

The fitted quantities are tabulated in Table 2. 

We use the global fit in Eq. (4.1) and depict in Fig. 1 and Fig. 2 the polarized proton 
structure function for some values of Q^ and compared them with experimental data. The 
comparison indicates that the global fit works well. 



Parameters 


valuse 


«o 


-0.0287 ±0.1463 


ai 


0.1253 ±0.0671 


«2 


-0.0092 ±0.0145 


a 


1.2624 ±0.0979 


b 


2.227 ±0.4745 


X^ (Goodness of fit) 


0.984 



Table 2: Global fit parameters obtained by fitting the Eq. (4.1) over the experimental data 
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Figure 1: Polarized proton structure function, gf for some values of Q^ extracted from global fit. 
Experimental data are from llsL p2]. The data cover the kinematic regions 0.0041 < x < 0.9 and 
O.lSGeV^ < Q^ < 20GeV^. 

4.2 Numerical results for the polarized gluon distribution 

In this section we will briefly describe the numerical results for polarized gluon distribution 
by using the analytic solution, Eq. (3.18). Finally we will compare our result for x6g{x, Q'^) 
with those from ^^C7'08, DSSV'08 and LSS'06 global fits pi. 



im- 



In this calculation, we use the LO approximation of as{Q^) which is defined in [|^ 



a.(Q') 



4tt 



/3i 



/3oln(Q2/A2)' 
2 



11 



-n 



/' 



(4.3) 
(4.4) 



with rif = 5 and A5 = 146 MeV for Q > 4.5 GeV, n/ = 4 and A4 = 192 MeV for 



xg, 



xgr 



0.04 



0.02 



COMPASS 2009 
HERMES 2006 
SLAC-E143 
globalf it 
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• 
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0.07 
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Figure 2: Polarized proton structure function, xg\ at Q^ — SGeV^ and Q^ — 5GeV^ extracted from 
global fit. Experimental data are from @ ^ |20[ |^, ||, ||l . 



1.3 GeV < Q < 4.5 GeV, and n/ = 3 and A3 = 221 MeV for Q < 1.3 GeV. These values 



have been used in CTEQ5L [26| and also used to extract the analytic unpolarized gluon 

distribution in [pl|]. 

Now by following these two steps we calculate the polarized gluon distribution: 

• Calculating 55(x,g2) by using the Eq. (2.2): gg(,x,Q^) = {^Y,e'^)~^Gf{x,Q'^). The 
convolution integrate in Eq. (2.2) with plus prescription, ()_(_, can be easily calculated 



using 1 17 1 



/ -fi-hg{y) = I -f{-My) - -g{x)]-g{x) f dyf{y) (4.5) 
Jx y y Jx y y y Jo 

• By using the Eq. (3.18), we can extract numerically the polarized gluon distribution 
inside the proton . 

Our result for polarized gluon distribution inside the proton is shown in Fig. 3 (left) 
and compared with some global fits. 

4.3 How to test the accuracy of the procedure ? 

An independent method of checking the numerical accuracy of the entire procedure for 
extracting the polarized gluon distribution is to go back to the original DGLAP equation 
from which we started, Eq. (2.3), i.e., 



fi f]^ ^ 

gg{x,Q^) = x -G{z,Q^)5K,g{-), 

Jx ^ ^ 



and numerically integrate its right hand side, which depends on our numerical solution . 
We then compare it with the gg{x,Q'^), which is independently known and arising from 
our global fit for the polarized proton structure function, g'f(x, Q^), based on its relation to 
GiixjQ"^) in Eq. (2.2). In Fig. 3 (right) we plot the results of these different methods. Test 
of accuracy shows that our solution for polarized gluon distribution is correct (Eq. (3.18)). 
Of course, we should note that for using this solution for polarized gluon distribution, 
we need to do a global fit with experimental data for g\{x,Q'^). We do not have enough 
experimental data for small x (x < 10~^) also for large x (x > 0.7), therefore we need to 
have (^{x^ Q"^) at these regions to be sure about the shape of polarized gluon distribution 
and future experiments should be focused on measuring the polarized structure function 
in these two regions. 
The gluon contribution to the spin of the proton in our calculation is 



A5(Q2 = bGeV'^) = I 5g{x, Q^ = hGeV'^)dx = 0.43 

JO.OOl 



(4.6) 



This value indicates that the gluon contribution to the proton spin is considerable and 



compatible with other phenomenological models |23, ^, 25] 
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Figure 3: Left: Polarized gluon distribution function, xSg{x,Q^) at Q^ = bGeV^ and comparison with 
other global fits. Right: Test of accuracy. 



5. Conclusion 

In this paper we tried to find an analytic solution for the DGLAP equation for polarized 
structure function and find the polarized gluon distribution inside the proton. Our solution 
is model independent and it is free of any parameters, initial input densities for solving the 
DGLAP equations and also free of the initial scale of Qq. It is only dependent on finding 
an analytic function for gf (x, Q^) by using the experimental data. Our results suggest that 
the precise measurement of experimental values for g\{x,Q'^) over large kinematic ranges 
of X and Q^ can predict directly the polarized gluon distribution more accurately. 
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A. The DGLAP equation for the polarized proton structure function 

In this part, we want to prove the Eq. (2.1). 

(A.l) 
We begin with the DGLAP equations for the polarized parton distribution functions 
for quark and ani-quark sectors 



%^''^ = Si/" 'i'i(^''^''<« (l) + /' T*»(-«^)^^« (1)1. (A.3) 

where t = LnQ'^ and we used the following properties for the polarized splitting func- 
tions |17|: 



5Kgg = 5Kqq, (A.4) 

5Kqg = 5Kqg, (A.5) 

After summation of the Eq. (A.2) and Eq. (A.3) and multiply both side by ^ X]i=i ^f^^ 
we have 

dg^{x,Q'^) as [^ dz ^2Nr7^ {^\ , "^ v^ 2 [^ dz r,2\xi^ P^ 

(A.6) 
In the Eq. (A.6) we have G{z, Q^) = z6g{z, Q^) and g\{x, Q"^) is the LO polarized spin 
structure function for proton 



Thus we achieved to our desired result in Eq. (2.1): 



(A., 
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